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possibility  of  developing  new  processing  strategies  for  high 
temperature  intermetaliic  compounds  is  being  investigated.  In  particular 
rapid  solidification  followed  by  controlled  heat  treatment  may  provide  new 
and  unusual  micro structures  of  multiphase  materials.  This  report  describes 
research  performed  at  NIST  to  develop  predictive  models  for  solubility 
extension  and  metastable  phase  formation  of  intermetaliic  compounds  and 
research  to  improve  the  phase  diagram  modeling  of  systems  involving  ordered 
phases,  t 

Theory  has  been  developed  to  predict  the  formation  of  disordered  form  of 
intermetaliic  compounds  during  rapid  solidification.  Using  a  modification  of 
the  Aziz  solute  trapping  theory,  solidification  velocities  required  to  form, 
for  example,  a  BCC  phase  from  the  melt  when  a  B2  phase  is  the  stable  phase 
have  been  developed.  In  the  present  reporting  period,  this  theory  has  been 
extended  to  the  case  where  a  FCC  phase  replaces  a  Ll0  phase. 

Experimental  research  has  continued  on  the  intermetaliic  system,  NiAl- 
NiTi,  using  laser  surface  melting  and  examination  by  transmission  electron 
microscopy.  In  the  NiAl-NiTi  system,  the  equilibrium  intermediate  Heusler 
phase,  Ni2AlTi,  was  suppressed  by  rapid  solidification  and  extension  ot  the 
composition  range  of  the  NiAl  and  NiTi  phases  was  observed  by  melt  spinning. 
To  quantify  the  solidification  rates  required  for  this  process,  samples  have 
been  prepared  by  pulsed  laser  melting  where  the  solidification  rate  can  be 
varied  between  1-10  m/s.  Preliminary  experiments  show  the  formation  of  the 
B2  phase  at  the  NizAlTi  composition,  which  under  equilibrium  processing  forms 
the  L2i  structure. 

Experimental  work  on  ternary  alloys  surrounding  the  composition  TizNbAl 
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to  determine  the  phases  present  in  arc  melted  and  heat  treated  samples  has 
continued.  TEM  investigations  are  required  to  sort  out  the  complex  ordering 
reactions  in  this  important  alloy  system.  A  broad  equilibrium  composition 
range  for  the  B2  phase  has  been  identified  at  1100°C  and  complex  low 
temperature  (~700°C)  transitions  are  being  examined. 

An  evaluation  of  existing  data  and  a  thermodynamic  calculation  of  the 
Ti-Al-Ta  phase  diagram  has  been  initiated  using  the  THERMOCALC  code.  This 
activity  was  requested  at  the  March  1988  review  of  the  DARPA/ONR  Program  on 
"Development  of  High  Temperature  Intermetallics  for  Structural  Aerospace 
Applications"  at  Pratt  &  Whitney,  Florida. 

I .  Introduction 

The  development  of  high  temperature  materials  is  closely  related  to  the 
formulation  of  processing  strategies  for  chemically  ordered  phases.  Most 
interraetallic  compounds  including  aluminides,  carbides,  and  silicides  as  well 
as  high  temperature  ceramic  phases  are  ordered.  However,  optimum  mechanical 
properties  are  likely  to  come  from  intimate  dispersions  of  several  phases, 
some  of  which  are  ordered.  These  dispersions  can  be  produced  by  a  phase 
transformation  sequence  involving  both  ordering  and  phase  separation  , 
beginning  with  a  solid  phase  of  a  carefully  selected  unstable  composition 
made  by  rapid  solidification. 

Recently,  significant  advances  have  occurred  in  the  utilization  and 
understanding  of  rapid  solidification  processing  of  alloys.  Factors  which 
promote  refined  segregation,  solubility  extension  and  metastable  phase 
formation  have  been  identified.  However,  much  of  this  research  has  been 
focused  towards  disordered  crystalline  phases:  i.e.,  terminal  solid 
solutions,  not  ordered  intermetal lie  compounds. 
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At  the  same  time,  significant  advances  have  been  realized  in  the 
thermodynamics  and  kinetics  of  order-disorder  transitions.  The  distinction 
between  first  and  higher  order  transitions  has  been  clarified,  the  kinetics 
of  ordering  reactions  and  the  structure  and  mobility  of  APB's  have  been 
determined  and  reactions  that  involve  fine  scale  ordering  and  compositional 
separation  have  been  studied. 

This  research  attempts  to  combine  the  advances  in  these  areas  to  develop 
new  processing  strategies  for  high  temperature  ordered  multiphase  materials. 

In  section  II  of  this  report  we  describe  theoretical  and  experimental 
research  focused  on  determining  the  possibility  of  extending  the  solubility 
range  of  ordered  phases  by  rapid  solidification.  Subsequent  heat  treatment 
of  these  metastable  alloys  can  form  stable  high  temperature  multiphase 
mixtures.  This  research  also  includes  an  examination  of  the  state  of 
nonequilibrium  order  of  rapidly  quenched  intermetallic  compounds. 

Section  III  of  this  report  summarizes  research  at  phase  identification 
in  composition  surrounding  TizNbAl. 

Section  IV  of  this  report  describes  phase  diagram  activities  involving 
ordered  phases.  A  preliminary  calculation  of  the  ternary  diagram  Ti-:Ta-Al 
has  been  produced. 

II .  Solubility  Extension  and  Disorderinz  of  Intermetallic  Compounds  by* Rapid 

Solidification 

(a)  Theory  of  Solute  &  Disorder  Trapping  bv  Rapid  Solidification- 

A  model  has  been  developed  to  predict  the  long  range  order  parameter  and 
composition  of  a  chemically  ordered  phase  as  a  function  of  interface 
velocity.  The  details  of  this  model  were  described  in  the  previous  semi¬ 
annual  report.  The  model  is  an  extension  of  the  solute  trapping  model  of 
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Aziz  and  relies  on  an  analysis  of  the  interdiffusion  across  the  liquid  solid 
interface  between  a  liquid  phase  and  a  solid  composed  of  two  sublattices. 

The  model  predicts  the  transition  from  solidification  of  a  solid  phase  with 
equilibrium  long  range  order  parameter  and  with  equilibrium  partitioning  of 
composition  at  low  interface  velocity  to  solidification  of  a  disordered 
crystalline  solid  with  the  same  composition  as  the  liquid  at  high 
solidification  velocity. 

A  broad  range  of  experimental  results  on  rapid  solidification  of 
intermetallic  compounds  suggests  that  solute  trapping  and  disorder  trapping 
do  not  occur  for  line  compounds  under  melt  spinning  conditions.  Trapping 
seems  only  to  occur  for  ordered  phases  in  which  the  ordering  is  not  so 
strong;  i.e.,  where  the  equilibrium  composition  range  is  large  (-10%).  The 
general  theory  was  explored  and  a  simple  result  related  to  this  point  has 
been  obtained.  The  critical  velocity  for  the  trapping  of  disorder,  Vc ;  i.e., 
the  velocity  above  which  only  the  disordered  variant  of  the  intermetallic 
compound  can  form  can  be  estimated  from 

TC 

V  =  V  (—  -  1) 

C  D'T  ' 
m 

where  VD  is  the  ratio  of  the  interface  diffusivity  to  the  jump  distance  (-102 
cm/s),  Tc  is  the  critical  temperature  for  the  order-disorder  transition,  and 
Tm  is  the  melting  point  of  the  compound.  When  the  compound  is  ordered  up  to 
its  melting  point  Tc/Tm  >  1.  One  can  see  that  the  higher  Tc ,  which  relates 
directly  to  the  strength  of  the  ordering,  the  higher  the  velocity  required  to 
trap  disorder.  For  line  compounds  Tc  »  Tm  and  thus  Vc  may  be  unattainable 
in  normal  rapid  solidification  methods. 

The  results  of  this  theory  have  been  calculated  for  solids  which  have 
first  and  second  order  ordering  transitions  during  the  current  reporting 
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period.  Discontinuities  in  both  order  parameter  and  solid  composition  occur 
when  the  growth  velocity  is  increased  when  the  equilibrium  ordering  reaction 
is  first  order.  These  results  are  described  in  detail  in  the  preprint 
enclosed  in  the  appendix. 

(b)  NiAl-NiTi  Alloys 

In  previous  reports  experimental  research  on  the  microstructure  of  melt 
spun  alloys  in  the  NiAl-NiTi  system  was  described.  Eutectic  alloys  of  NiAl- 
Ni2TiAl  and  Ni2TiAl-NiTi  solidified  as  a  single  B2  phase.  However,  results 
on  alloys  with  composition  Ni2TiAl  were  not  clear.  The  APB  size  was 
sufficiently  coarse  that  the  non -equilibrium  solidification  product  could  not 
be  determined  with  certainty.  In  collaboration  with  Professor  Mike  Aziz  of 
Harvard  we  have  performed  pico-second  pulse  melting.  This  technique  produces 
solidification  rates  of  1-5  m/s.  Initial  results  show  that  extremely  fine 
L2:  domains  are  produced  by  this  technique  suggesting  that  the  solidification 
process  produced  the  B2  phase  with  subsequent  solid  state  ordering  to  the  L2j 
structure.  Future  experiments  with  this  method  will  permit  testing  of  the 
theory  described  in  the  previous  section. 

III.  B2  and  Related  Phases  in  the  Ti-Al-Nb  System 

Titanium  aluminides  (Ti3Al  and  TiAl)  with  -10  at%  Nb  additions  have 
received  considerable  attention  as  potential  low  density,  high  strength’  and 
creep  resistant  materials.  However,  the  phase  equilibria  in  this  ternary 
system  is  poorly  understood.  Because  the  presence  of  BCC-based  phases  in 
these  alloys  seems  to  play  an  important  role  in  the  deformation,  in-depth 
studies  of  the  BCC-based  phase  fields  are  being  conducted.  A  series  of 
alloys  surrounding  the  composition  Ti2AlNb  are  being  studied  by  TEM  to 
determine  the  structure  of  phases  present  in  arc  cast,  melt  spun  and  heat 
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treated  samples.  When  complete,  these  results,  combined  with  those  from  the 
University  of  Wisconsin  by  J.  Perepezko  and  Y.  A.  Chang  under  DARPA 
sponsorship,  will  be  combined  with  ternary  phase  diagram  calculations 
performed  under  this  contract  to  determine  the  ternary  phase  diagrams  Ti-Nb- 
Al. 

(a)  Literature  Survey 

Experimental  work  focused  on  ternary  equilibria  in  the  BCC  field  of  the 
Ti-Al-Nb  system  is  sparce.  Strychor,  Williams  and  Soffa  (Met.  Trans  19A 
(1988)  225)  have  shown  that  alloys  with  -25  at%  Al  and  between  7  and  20  at% 

Nb  which  are  quenched  from  1250°C  are  single  phase  B2  with  APB's.  This 
indicates  that  the  BCC  to  B2  transition  for  these  compositions  occurs  below 
1250°C.  Quenched  alloys  with  0-7  at%  Nb  are  HCP  a'  martensite.  The  B2 
alloys  also  contain  a  tweed  structure  of  ordered  o>  phase.  During  aging  at 
400°C  for  100  hours,  the  B2  phase  was  partially  consumed  by  an  ordered  w 
phase.  Evidence  for  additional  ordering  of  w  into  a  B82  (ZrzAl)  phase  is 
also  presented.  No  D019  phase  was  observed  to  form  for  alloys  with  greater 
than  7  at%  Nb. 

Banerjee,  Gogia,  Nandi  and  Joshi  (Acta  Met.  36  (1988)  871)  have  studied 
the  composition  Ti-25  at%  Al-12.5  at%  Al.  Alloys  held  at  llOO'C  fbr  3  hours 
and  quenched  contain  a2(D019)  and  B2  phases.  The  B2  phase  contains  no  APB's. 
This  indicates  the  equilibria  at  1100°C  is  D0ig  +  B2,  not  D019  +  BCC. 

Furnace  cooling  of  this  structure  yields  microstructures  that  suggest  that 
both  of  the  phases  transform  partially  to  an  orthorhombic  0(Cd3Er)  phase  with 
symmetry  Cmcm.  Evidence  is  presented  that  a2  •*  a z  +  0  by  simultaneous 


ordering  and  spinodal  decomposition  and  that  the  B2  phase  also  transforms  by 
nucleation  and  growth  to  the  orthorhombic  phase.  They  also  infer  that  the 


orthorhombic  phase  has  a  preferred  stoichiometry  of  Ti2NbAl. 

Banerjee,  Nandy  and  Gogia  (Scripta  Met.  21.  (1987)  597)  have  examined  the 
site  occupancy  of  the  B2  phase  which  occurs  in  a  two-phase  mixture  with  the 
a2  phase  in  heat  treated  samples  (no  temperature  given)  of  Ti-25.6  at%  Al- 
10.1  at*  Nb.  Using  channeling  enhanced  microanalysis  they  determined  that 
one  sublattice  was  mostly  Ti  while  the  other  contained  the  excess  Ti  and  the 
Al  and  Nb  atoms.  Using  site  fractions  x) ,  they  give 

i  xTi  XA1  XNb  i  _  -96  0  .04 

I  !  ' 

I  x?i  xai  x£b  j  -28  -48  24 

(b)  Experimental  Procedure 

Alloys  shown  in  Table  1  were  prepared  by  arc  melting.  In  excess  of  ten 
remelts  were  necessary  to  homogenize  the  buttons. 

Table  1 .  Ti-Al-Nb  Alloys  Under  Investigation 


Ingot  * 

at*  Ti 

at*  Al 

at*  Nb 

ppm  0 

ppm  N 

ppm  H 

30 

50  (51.1) 

25  (23.9) 

25  (25.0) 

500 

130 

22 

01 

50  (50.6) 

12.5  (12.2) 

37.5  (37.2) 

500 

320 

34 

02 

37.5  (38.1) 

25  (24.3) 

37.5  (37.6) 

290 

30 

7 

03 

62.5  (63.5) 

25  (24.2) 

12.5  (12.3) 

520 

90 

12 

04 

50  (51.4) 

37.5  (36.2) 

12.5  (12.4) 

630 

50 

7 

Cast  structures  were  homogenized  with  a  3  hour  heat  treatment  at. 1400°C 
in  a  vacuum-tight  furnace  under  2/3  atm  gettered  Argon.  It  was  found 
necessary  to  rest  samples  on  a  Y203  coated  Al203  substrate  to  prevent 
reaction  with  A1203 .  Cooling  was  performed  by  lowering  the  samples  out  of 
the  hot  zone  of  the  furnace  into  a  lower  chamber.  The  cooling  rate  was 
estimated  by  visual  observation  and  by  measurements  of  the  furnace  vendor  to 
be  -400  K/min.  Heat  treatments  at  1100°C  were  performed  after  the  1400°C 
treatment  by  lowering  the  furnace  temperature.  Equilibration  at  1100°C 
occurred  in  -5  min.  Heat  treatments  at  700°C  were  performed  by  sealing 
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individual  samples  from  the  1400°C  treatment  in  quartz  tubes.  Samples  were  j 

wrapped  in  Ta  foil  and  tubes  were  backfilled  with  2/3  atm  He.  ' 

Chemical  analysis  was  performed  on  the  samples  homogenized  at  1400 °C. 

Because  these  samples  were  nominally  single  phase  when  viewed  by  optical  or  i 

SEM  metallography,  microprobe  analysis  with  elemental  standards  yields 
reliable  analysis  of  the  bulk  alloy  compositions.  The  values  obtained  are 
given  in  parentheses  in  Table  1.  Oxygen,  nitrogen  and  hydrogen  analysis  was 
performed  on  these  samples  by  a  commercial  vendor  and  is  also  given  in  Table 
1. 

(c)  Results 

The  results  of  microstructural  investigation  of  phase  identification  to 
date  are  summarized  in  Table  2.  Considerable  solid  state  transformation 
occurs  during  cooling  in  some  samples.  Actual  phases  present  are  listed 
along  with  the  inferred  equilibrium  at  temperature  and  the  decomposition 
path.  The  results  combine  observations  made  by  optical  metallography,  x-ray 
diffraction  and  electron  microscopy. 

Microstructure  of  Arc  Melted  Buttons  -  Optical  metallography  of  the  as- 
cast  buttons  revealed  dendritic  but  single  phase  structures  in  all  samples 
except  #03.  Microsegregation  was  small  (<1%)  but  the  cores  of  the  dendrites 
were  rich  in  Nb  as  judged  from  backscatter  contrast.  Observation  of  grain 
boundaries  cutting  across  dendrite  arms  suggests  significant  grain  coarsening 
(up  to  300  ^im)  during  solid  state  cooling.  In  sample  #03,  a  fine  martensitic 
plate  structure  was  visible.  X-ray  diffraction  indicated  the  structure  to  be 
D019 (Ti3Al) . 

Microstructure  of  Samples  Heat  Treated  at  1400°C  -  Optical  metallography 
of  the  samples  cooled  from  1400°C  reveals  the  absence  of  dendritic 
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microsegregation  and  very  coarse  grains  (-2  mm) .  The  microstructures  appear 
essentially  featureless  with  the  exception  of  #03  which  is  martensitic.  TEM 
examination  of  the  martensitic  sample  indicates  complete  transformation  from 
the  high  temperature  phase  which  existed  at  1400°C  to  a  D019  martensite.  The 
other  four  compositions  consist  primarily  of  a  B2  or  related  phase.  Although 
no  APB's  were  observed,  previous  research  by  Banerjee  et  a 1.  and  Strychor  et 
al.  suggest  a  BCC  -»  B2  ordering  occurs  -1200°C.  We  thus  infer  that  the 
equilibrium  phase  at  1400°C  is  BCC  for  all  five  alloys.  However,  the 
possibility  exists  that  for  some  compositions  the  B2  is  stable  at  1400°C. 

High  temperature  x-ray  diffraction  will  be  used  to  answer  this  point. 

Besides  the  B2  ordering  which  occurs  during  cooling,  sampLe  #04 
completely  transforms  to  ordered  o  (or  possible  B8Z ) .  Sample  #30  contains 
small  plates  of  twinned  D019. 

Microstructure  of  Samples  Heat  Treated  at  1100°C  fcr  Four  Days  -  The 
microstructure  of  samples  #03,  #30,  and  #01  cooled  from  1100°C  are  identical 
to  those  cooled  from  1400°C.  Oddly  the  result  for  #03,  that  the  phase 
present  at  1100°C  after  3  hours  is  B2 ,  directly  contradicts  the  result  of 
Banerjee  et  al.  who  found  B2  +  D0ig  under  these  conditions.  This  may  be 
related  to  differences  in  cooling  rate. 

Sample  #02  contained  B2  and  a  phases  while  sample  #04  contained  B2  + 
D0lg.  The  B2  in  this  alloy  #04  was  completely  transformed  to  ordered  u. 

(d)  Future  Research 

Future  work  will  involve  an  analysis  of  the  phases  present  at  700°C 
after  a  one  month  annealing.  The  calculated  ternary  phase  diagram  Ti-Al-Nb 
which  was  presented  in  the  last  report  will  be  modified  to  be  consistent  with 
this  new  data  and  other  data  generated  by  J .  H.  Perepezko  and  Y.  A.  Chang  at 
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the  University  of  Wisconsin,  Madison,  under  DARPA  sponsorship. 

IV .  Phase  Diagram  Modeling 

We  have  undertaken  a  preliminary  calculation  of  the  Ti-Al-Ta  ternary 
system  in  response  to  requests  by  Ralph  Hecht,  Pratt-Whitney ,  Florida,  for 
use  in  their  DARPA  Program  "High  Temperature  Metallics  for  Structural 
Aerospace  Applications." 

A  preliminary  calculation  of  the  Ti-Al-Ta  ternary  phase  diagram  was 
undertaken  using  the  THERMOCALC  DATABANK  system  described  by  B.  Sundman,  B. 
Jansson  and  J.  0.  Anderson  (CALPHAD  9,  2  (1985)  pp.  153-190).  The  stability 
ranges  were  adjusted  to  match  the  isothermal  section  at  1100°C  by  Sridharan 
and  Nowotny  (Z.  Metallkunde  1J±  (1983)  468). 

In  the  calculations  of  the  binary  systems  the  liquid  and  elemental 
phases  were  described  as  sub-regular  solutions.  The  intermetallic  compounds 
in  the  Ta-Al  and  Ti-Al  were  described  with  the  Wagner -Schottky  model,  where 
intermetallic  compounds  are  considered  to  consist  of  different  sublattices 
allowing  substitutional  solutions  on  each  of  these  sublattices.  Since  the 
Wagner -Schottky  model  is  mathematically  a  special  case  of  the  sublattice 
model,  the  parameters  of  these  intermetallic  compounds  can  be  easily  , 
transformed  to  the  sublattice  model  description  of  the  THERMOCALC  DATABANK 
system.  For  the  present  only  the  a  and  TaAl3  intermetallics  were  modeled  in 
the  Ta-Al  binary. 

For  the  extrapolation  of  ternary  Ti -Ta-Al  system  all  phases  must  be 
allowed  to  have  ternary  ranges  of  homogeneity,  which  means  for  the  Wagner - 
Schottky  phase  that  they  have  to  be  modeled  as  metastable  phases  for  the  two 
other  binary  systems. 

As  a  first  approximation  the  Gibbs  energies  of  formation  of  the  Ta-Al 


12 


compounds  in  the  Ti-Al  system  and  the  Ti-Al  compounds  in  the  Ta-Al  system 
were  assumed  to  be  the  same  as  in  the  system  where  these  compounds  are 
stable,  for  the  Ta-Ti  system  the  Gibbs  energies  of  formation  were  assumed  to 
be  the  same  as  for  these  compounds  consisting  only  of  Ta  or  Ti,  respectively. 

In  order  to  prevent  these  metastable  phases  from  becoming  stable  in  the 
binary  systems  and  to  adjust  the  calculated  phase  diagram  to  experimental 
results  at  1100°C,  the  Gibbs  free  energies  of  formation  of  these  metastable 
phases  were  adjusted  by  assuming  that  the  entropies  of  formation  of  the 
metastable  phases  have  the  same  magnitude  as  the  stable  phases  and  the 

enthalpies  of  formation  of  the  metastable  phases  are  more  positive  than  those 

of  the  stable  phases.  The  ternary  ranges  of  homogeneity  at  1100°C  calculated 

with  these  parameters  were  acceptable  except  for  the  hep  and  Ll0  phase,  where 

they  were  too  small.  This  was  adjusted  by  giving  the  Ta-Al  interaction  term 
of  the  hep  phase  the  most  negative  value  possible,  without  the  hep-phase 
becoming  stable  in  the  Ta-Al  system.  For  the  Ll0 -phase  ternary  interaction 
terms  were  introduced  between  Al  and  Ta  on  the  first  sublattice  and  with  Ti 
on  the  second  one  and  between  Al  and  Ti  on  the  first  sublattice  and  Ta  on  the 
second  one. 

The  calculated  binary  diagrams  used  for  the  Ti-Ta-Al  ternary  ^calculation 
are  shown  in  Figures  1,  2,  and  3.  The  correct  Ti-Al  and  Ta-Al  diagrams’  are 
still  subject  to  debate  in  some  regions  and  can  be  modified  at  a  later  time. 

The  results  obtained  from  the  calculations  with  these  parameters  are 
shown  in  Figures  A,  5,  and  6.  Ternary  phases  were  omitted  at  this  initial 
stage.  Figures  A  and  5  show  isothermal  sections  at  1100,  and  1A00°C,  while 
Figure  6  shows  the  liquidus  projection. 

Further  work  on  phase  diagram  modeling  in  this  system  must  await  a  more 
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precise  definition  of  the  measured  phase  boundaries. 


V.  Fiscal  Status 

(1)  Amount  currently  provided  for  contract  program 

$150,000  for  period  February  9,  1987  to  September  30,  1987. 

$200,000  for  period  October  1,  1987  to  September  30,  1988. 

(b)  Expenditure  and  commitment  during  April  1,  1988  to  September  30,  1988 

$100,000 

(c)  Estimated  funds  required  to  complete  this  work 


October  to  December  1988  $50,000 
January  to  March  1989  $50,000 
April  to  June  1989  $50,000 
July  to  September  1989  $50,000 


October  to  December  1989 


$50,000 


Figure  Captions 


Figure  1.  Calculated  Ti-Al  phase  diagram  used  for  calculation  of  Ti-Al-Ta 
ternary  diagram. 

Figure  2.  Calculated  Ti-Ta  phase  diagram  used  for  calculation  of  Ti-Al-Ta 
ternary  diagram. 

Figure  3.  Calculated  Ta-Al  phase  diagram  used  for  calculation  of  Ti-Al-Ta 
ternary  diagram. 

Figure  4.  Calculated  isothermal  section  of  the  Ti-Al-Ta  ternary  phase 
diagram  at  100°C  omitting  ternary  phases. 

Figure  5.  Calculated  isothermal  section  of  the  Ti-Al-Ta  ternary  phase 
diagram  at  1400°C  omitting  ternary  phases. 

Figure  6.  Calculated  liquidus  and  solidus  invariants  of  the  Ti-Al-Ta  ternary 
phase  diagram  omitting  ternary  phases  (liquidus  invariants  have  double 
arrows) . 
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THERMO-CALC  (87.12.24:12.17)  :Ti-Al 


MOLE-FRACTION  AL 


THERMO-CALC  (88.10.  6:17.  7)  :Ti-Ta 


MOLE-FRACTION  TA 


THERMO-CALC  (83.10.12:15.30) 


MOLE-FRACTION  TA 


THERMO-CALC  (88.10.12:1 


MOLE-FRACTI 


THERMO-CALC  (88.10.13 


u 


MOLE-FRA 


APPENDIX 


A  Theory  for  the  Trapping  of  Disorder  and  Solute  in  Intermetallic 
Phases  bv  Rapid  Solidification 

William  J.  Boettinger 
Metallurgy  Division 

National  Institute  for  Standards  and  Technology 
(formerly  National  Bureau  of  Standards) 

Gaithersburg,  MD  20899 

and 

Michael  J.  Aziz 
Division  of  Applied  Sciences 
Harvard  University 
Cambridge,  MA  02138 

A  theory  is  developed  to  predict  the  long  range  order  parameter,  composition  and 
interfacial  temperature  of  a  chemically  ordered  phase  as  a  function  of  interface  velocity 
during  rapid  crystal  growth.  The  theory  extends  the  solute  trapping  theory  of  Aziz  to  a 
solid  phase  consisting  of  two  sublattices.  The  engulfment  of  atoms  randomly  on  the  two 
sublattices  by  the  rapidly  moving  interface  is  balanced  against  the  interdiffusion  across 
the  liquid-solid  interface  which  attempts  to  restore  local  equilibrium.  With  increasing 
interface  velocity  the  theory  predicts  a  transition  from  the  solidification  of  a  phase  with 
equilibrium  long  range  order  parameter  and  with  equilibrium  solute  partitioning  to  the 
solidification  of  a  disordered  crystalline  phase  with  the  same  composition  as  the  liquid. 
Predictions  for  various  free  energy  functions  for  the  solid  phase  suggest  that  the  decrease 
of  order  parameter  with  increasing  interface  velocity  may  be  continuous  or  discontinuous 
and  that  transitions  to  solute  trapping  and  to  disorder  trapping  can  occur  at  different 
growth  rates. 

1.  Introduction 

Over  the  years,  considerable  attention  has  been  paid  to  the  kinetics  of  motion  of  the 
interface  separating  two  phases,  especially  to  the  kinetics  of  solidification  of  a  melt. 


Recent  interest  in  the  rapid  solidification  processing  of  alloys  and  other  materials  using 
transformations  far  from  equilibrium  has  renewed  interest  in  these  kinetics.  Theories  of 
nonequilibrium  incorporation  of  solute  into  rapidly  growing  crystals  (1-9)  are  currently 
the  focus  of  experimental  tests  (10,11).  An  interesting  additional  possibility  occurs  for 
intermetallic  compounds  or  other  crystals  with  long  range  chemical  order.  Rapid  growth 
can  also  cause  the  formation  of  crystals  with  nonequilibrium  long  range  order.  Because 
ordering  in  the  solid  state  during  post-solidification  cooling  can  mask  events  at  the 
liquid-solid  interface  quantitative  experiments  to  measure  the  state  of  long-range  order  of 
the  solidifying  material  are  difficult  to  design.  However,  the  observation  of  solidified 
phases  with  high  densities  of  antiphase  domains  when  the  phase  is  normally  ordered  to  the 
melting  point  is  a  clear  indication  that  disorder  has  been  trapped  by  the  process  of  rapid 
solidification  (12-15). 

Kinetic  theories  for  trapping  generally  fall  into  two  categories:  diffuse  interface  and 
chemical  reaction  rate  theories.  Diffuse  interface  theories  (16-17)  solve  the  Cahn-Hilliard 
Equation  in  a  moving  frame  and  predict  in  a  formal  way  that  certain  kinds  of  disorder  can 
be  induced  into  a  growing  phase  by  rapid  interface  motion.  Chernov  has  developed  a 
chemical  reaction-rate  theory  that  has  been  applied  to  both  solute  trapping  during  rapid 
solidification  of  disordered  solid  solutions  and  disorder  trapping  during  solidification  of 
stoichiometric  compounds  (18).  Analytic  and  Monte-Carlo  models  of  solute  trapping 
based  on  very  similar  principles  to  Chernov's  have  been  developed  by  others  (19-22). 
Probabilities  are  assigned  for  each  species  to  hop  into  or  out  of  the  crystal  based  on  the 
temperature  and  the  chemical  (and  structural,  e.g.,  location  at  a  kink)  environment;  the 
net  hopping  rates  are  summed  to  yield  the  growth  velocity.  In  the  language  of  the 
thermodynamics  of  irreversible  processes,  the  flux  of  each  species  across  the  interface  is  a 
function  of  the  conjugate  driving  force  (chemical  potential  difference  across  the  interface) 


2 


for  that  species  only,  being  independent  of  the  driving  force  for  the  other  species. 
According  to  these  models,  for  the  impurity  to  be  incorporated  into  a  high-energy  site 
(the  wrong  sublattice  in  disorder  trapping;  any  lattice  site  in  solute  trapping),  it  must 
actively  hop  into  it. 

The  main  innovation  in  the  approach  of  Aziz  (1)  is  that  if  the  impurity  does  nothing 
it  may  end  up  on  a  high-energy  site  by  virtue  of  the  formation  by  its  neighbors  of  a 
regular  lattice  around  it.  “Hence  to  avoid  incorporation  onto  a  high-energy  lattice  site,  an 
atom  must  diffuse  away.  (For  growth  of  compounds  some  of  the  atoms,  of  course,  will  by 
chance  be  incorporated  onto  the  right  sublattice  even  if  they  do  no  hopping.)  Since  the 
maximum  speed  of  diffusion  can  be  rather  slow  compared  to  the  speed  with  which 
crystal-melt  interfaces  have  lx*en  observed  to  move  (23-25),  the  atom  may  be  trapped  on 
a  high-energy  site  by  a  rapidly  moving  interface.  This  implies  a  strong  cross-coupling 
between  the  fluxes  and  driving  forces  of  the  individual  species  (26).  These  assumptions 
have  been  shown  experimentally  to  be  much  more  reasonable,  at  least  for  solute  trapping 
during  rapid  solidification  (10-11),  than  that  requiring  each  atom  to  actively  hop  into  its 
final  site  in  the  crystal. 

In  this  paper  we  develop  a  chemical  reaction  rate  theory  for  the  kinetics  of 
crystallization  of  an  ordered  phase  at  an  atomically  "rough"  interface;  i.e.,  there  are  no 
special  sites  such  as  ledges  or  kinks.  The  interface  temperature,  composition  and  order 
parameter  of  the  solid  will  be  described  as  functions  of  the  interface  velocity  and  liquid 
composition  at  the  interface.  It  is  the  first  theory  to  treat  the  simultaneous  trapping  of 
solute  and  disorder.  This  is  carried  out  by  modeling  solute  trapping  of  the  undesirable 

*Yltf  A  TA 

species  on  each  sublattice.  The  "Sima,  of  the  sublattice  compositions  yields  the  overall  solid 
composition;  their  difference  yields  the  order  parameter.  First,  ohe  conditions  for 
equilibrium  between  a  liquid  and  a  solid  phase  consisting  of  two  equivalent  sublattices  are 


developed  in  a  form  required  for  the  kinetic  theory.  Second,  the  kinetic  theory  is 
presented;  and  third,  numerical  results  for  the  growth  of  solid  phases  with  First  and  second 
order  order-disorder  transitions  are  presented. 

2.  Thermodynamics  of  the  Phases  at  Equilibrium 


We  assume  for  simplicity  that  the  solid  phase  consists  of  two  equivalent  sublattices 

T 

o  and  0  with  no  vacancies  and  that  a  total  free  energy  function  for  the  solid,  G  is  a 
function  of  nj,  the  number  of  moles  of  component  i  (i  =  A,B)  on  sublattice  j  (j  =  a,0). 
The  constraint  of  equal  number  of  lattice  sites  of  the  two  sublattices  requires  that  n «  + 
iin  =  One  can  formally  define  the  chemical  potential  of  component  i  on  the  j 

sublattice  as 


chi  ■) 
1 


[1] 


The  conditions  for  equilibrium  in  the  solid  phase  are  given  by 

a  a  J  3 
/*B  /«a  “  ;/B  "/V 


If  we  assume  that  the  constraint  on  the  number  of  moles  on  each  sublattice  also  holds  in 


the  solid  at  an  interface  with  a  liquid  phase  then  one  obtains  at  equilibrium 

2  K  +  ;ia)  =  /(A 

2  ^'B  +  =  f‘B  ^ 

where  and  pj^  are  the  chemical  potentials  of  the  liquid  phase.  As  a  consequence  of  [2] 

and  [3],  the  conditions 

a  a  _  L  L 
/iB  MA  -  /'B  /‘a 

0  0  11 

UB  flA  =  flB~V\  [41 

also  hold  at  equilibrium.  These  conditions  must  be  satisfied  at  the  zero  growth  velocity 
limit  of  any  kinetic  theory. 

s 

Flic  free  energy  of  the  solid  phase  is  usually  written  on  a  molar  basis,  G  ,  in. terms 


3 


Ot  Ck  3  0 

of  the  compositions  of  A  and  B  atoms  on  the  n  and  0  sublattices,  x^,  Xg.  x  ^ ,  xj^  given 
bv 


a  S  ,  1 

XA  =  XA  +  2  * 

(i  S  1 

xA  =  xA~2f? 

a  S  1 

xB  =  xB-2? 
S  ,  1 

=  XB  +  2  r?‘ 


K0 

*B 


The  long  range  order  parameter,  t]  is  given  by 

*A 
S 


7?  =  xa-xa  =  xb 


s  s 

and  the  solid  compositions,  x  . ,  Xg,  are  given  by 

n  i 


B 


[6] 


yS  1(  0  _L  v@\ 

XA  ~  2*XA  +  XA> 
y  S  _  1/  »  , 

XB  ~  2'XB  +  xb' 


where  clearly.  xa  +  xb  ~  i,  x^  -r  *g  =  »,  x^  -r  XB 

c 

composition,  the  molar  free  energy  Gm  is  given  by 


['] 

3  S  S 

g  =  1,  x^  +  Xg  =  1.  With  these  definitions  of 

[S] 


Gm  _  2^*AXA  +  ^BXB  +  ^AXA  +  ^BXB^ 


ni  inn  u  u  r\  r\  u  u 

C  /3 

If  is  treated  as  a  function  of  independent  variables  Xg  and  Xg 

dGS 

a  a  c,  m 
^B  ~  ^A  =  2  ITTot 


P  then 


P  p 
^A 


dx 

OG 


=  2 


B 

S 

m 


dx 
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[10] 


and 


DC  ^ 

1  /  a  ,  ,  (i\  S  a  m 

2  <"A  +  "A>  =  Gm  -  XB  -TS 

B 
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-  x 


dx 

OG 


p  dG 


[H] 
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2  +  f4)  “  Gm  +  (1  “  XB^  -  XB>  T~J-  i12l 

dxB  dxB 

c 

If  on  the  other  hand  independent  variables,  Xg  and  77,  are  used 
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a  ,  ft  _  m 

3  0  < 

''B  “  ''A  =  77 


1  t  »  .  0\  _  r^S  S 

2  +  ^  ~  Gm  XB 


q  dGS  0GS 
S  in  in 

o\B 

'V'l  S  ATI  i 

S\  ^m 


.  ft  p  n  l/VJ  t/VJi  _ 

2  +  ^b)  =  Gm  +  ^  _xii^  T~?~  ’  l12^ 

5xb 

It  can  now  be  seen  that  Eqn.  [2]  is  equivalent  to  the  ordinary  condition  for  equilibrium 

0C,H 

or(lor-  ~TT  =  °- 

3.  Thermodvnamic  Constraint  on  Solid  Composition  and  Order  Parameter  During 


Crystallization 

The  solidification  process  must  always  involve  a  decrease  in  free  energy  AG.  Baker 
and  Calm  (27)  have  described  the  domain  of  possible  solid  compositions  that  can  form  for 
various  liquid  compositions  at  a  given  temperature  by  requiring  that  AG  <  0  where 

AG  =  G®  -  +  ,&SB)  |13] 

which  is  Eqn.  [15]  of  Ref.  27  in  a  different  form.  Eqn.  [13]  is  the  basis  for  the  usual 
"tangent  to  curve  rule"  to  graphically  show  the  change  in  free  energy  for  a  given  phase 
change.  Here  Gm  is  evaluated  for  the  composition  of  the  solid  phase  and  and  pB  are 
evaluated  for  the  composition  of  the  liquid  phase.  Application  of  this  principle  to  the 
present  case  will  yield  domains  of  possible  order  parameter  and  solid  composition  that  can 
solidify  from  various  liquid  compositions  at  a  given  temperature.  Any  prediction  of  a 
kinetic  model  must  yield  values  of  t),  solid  composition,  temperature  for  each  liquid 
composition  which  will  give  AG  <  0. 
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Figure  l  shows  schematically  the  domain  of  possible  order  parameter  which  can 


form  at  various  temperatures  for  a  congruently  melting  compound  if  the  liquid 
composition  is  equal  to  the  solid  composition.  These  domains  as  well  as  the  detailed 
kinetic  results  described  later  in  this  paper  depend  on  whether  the  ordered  phase  differs 
from  the  disordered  phase  by  a  second  or  first  order  transition  depicted  in  (a)  and  (b), 
respectively.  At  the  top  are  shown  the  free  energies  of  the  solid  and  liquid  at  three 
temperatures.  At  the  equilibrium  melting  temperature  Tp  only  one  value  of  rj  is  possible. 
Below  Tj,  AG  <  0  defines  a  range  of  possible  rj  given  in  the  cross  hatched  region  in  the 
bottom  of  the  figure.  Note  that  only  when  the  T  <  in  (a)  or  T  <  T.^  in  (b)  can  a  solid 

S.F-  •! 

with  r)  =  0  form.  In  addition  to  the  range  of  tj  which  may  be  possible,  a  range  of  4+qwd 
composition  for  each  i]  ami-  solid  composition  is  also  described  by  the  condition  AG  <  0 
(not  shown). 

4.  A  Kinetic  Model  for  the  Interface 

4.1  The  Solidification  Velocity-Composition  Relation 

We  treat  here  the  case  of  steady  state  continuous  growth  of  an  atomically  rough 
planar  liquid-solid  interface.  Following  the  work  of  Aziz  (1)  for  random  crystalline  solids, 
the  atomistic  interdiffusion  flux  between  adjacent  liquid  and  solid  monolayers  must  have 
a  specific  value  related  to  the  growth  rate  and  the  compositions  of  the  monolayers  in  order 
to  maintain  the  steady  state.  For  the  growth  of  an  ordered  phase,  the  interdiffusion  flux 
between  the  liquid  and  each  sublattice  must  be  considered. 

To  treat  interdiffusion  between  the  liquid  and  solid  phase  during  crystallization  we 

s' 

use  redistribution  potentials  (/Zg  -  /z^)  as  defined  by  Aziz  and  Kaplan  (2)1;  viz., 


'Redistribution  potentials  can  also  be  written  in  terms  of  activity  coefficients  7^,  7^  and 


free  energies  of  the  pure  components  G^,  GR  as  (//R  -  fi^'  =  GR  -  G^  +  RT  fri  7g  - 

V' 

■  ‘  .  *  ■"  :i  * 

/-*  * 
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# 


(ab  -  /<A)  =  /'B  ~^A~KT  (n  XB  +  RT  61  XA'  1141 

I  L  ^ 

Redistribution  potentials  for  the  liquid  phase,  (/r^  -  flA)  >  and  the  a  and  0  sublattices  of 
the  solid  phase,  (//^  -  /<A)  and  ( ft ^  -  /*^)  arc  obtained  from  a  thermodynamic  model  of 
the  phases. 

We  consider  a  pair  of  reactions  from  initial  states  (of  B  in  an  o  or  0  site  in  the  solid 

adjacent  to  A  in  the  liquid)  over  a  barrier  to  final  states  (where  the  two  atoms  have 

exchanged  positions)  as  shown  in  Figure  2.  In  a  reference  frame  fixed  on  the  crystal 

lattice  the  forward  reactions  are  designated,  J^a  and  and  the  reverse  reactions, 

and  and  will  be  expressed  in  exchanges  per  unit  time  per  unit  area  of  a  (or  (i) 

sublattice.  We  assume  the  barrier  to  redistribution  remains  a  constant  height  above 

the  state  (initial  or  final)  with  the  higher  redistribution  potential.  Although  other 

assumptions  are  possible  for  the  definition  of  the  barrier  height  this  choice  is  consistent 

with  that  used  for  continuous  growth  of  pure  melts  (28).  Transitions  from  higher  to  lower 

V 

redistribution  potential  are  written  as  the  product  of  an  attempt  frequency  \  with  a 


% 

Boltzman  factor,  exp  {-  ypp}-  Transitions  from  a  lower  to  a  higher  redistribution 
potential  are  written  as 


/exp  hqd  + 1  (4  -  4)"  -  (4  -  4)1/rt> 

The  state  which  is  higher  or  lower  depends  on  the  values  of  order  parameter  and 
composition. 

Due  to  the  above  assumptions  regarding  the  barrier  heights  it  is  necessary  to 
introduce  parameters  Aj  (i  =  a,0)  with 

°if(4~4)  ~(4~4)  - 0 


1  if  (4  ~  4^  -  (4  <  0  • 


[15] 


RT  fa  7A. 
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Thus  the  interdiffusion  fluxes  are  given  by 

MQd  +  A“|(/.g  -  „%)' 


,  +  rv  _  f//A  an  L, 

JD  --jy-Xg(l  xg/  0XP 


,  L  L> 
(Mg  P  ) 


.V 


ITT 


JD°“TrxB^1“xB)exp 


[-(QD  +  (1  -  Aa)IO‘g  -  vtf  -  (/£  -  k)\\ 

RT - 


,+f]  _  iv\  3n  „L, 

'  D  -TrxB(l  xB)exp 


[Qd  +  A/jl(Mg  M^)  (Mg  m]() 
RT - 


i— . 3  _  f//A  L,.  /?\ 

'  D  ~  TTXB(1  “XB^  exp 


[Qd  +  (i  A^)|(/tg  fij^)  (hq  /i^) 


TFT 


[16] 


The  parameter  f  is  the  fraction  of  sites  at  the  interface  where  jumps  can  occur,  A  is  the 
jump  distance,  and  0  is  the  atomic  volume  (assumed  the  same  for  A  and  B).  This 
assumption  regarding  the  barrier  heights  is  identical  to  that  made  by  Aziz  and  Kaplan 
(2).  It  implies  that  the  exchanges  are  governed  by  the  properties  of  the  interface  and  not 
the  bulk  solid  or  liquid  phases.  One  also  could  have  introduced  separate  Qg  and  for 
the  exchanges  from  the  two  sublattices,  the  consequences  of  which  will  be  deferred  to  the 
discussion  section  of  this  paper. 

The  net  diffusive  fluxes,  Jg,  and  Jg,  are  then 

,o  _  ,+ar  T— a 
JD~JD  JD 


A.  T  ( 1  -  A  ) 

L/  n  v  n- 


=  D./(An)[xg(l  -  x^)Ka  a  -  x^(l  -  xg)Ka  a  ] 


[17] 


and  Jg  = 


r>i/(A!J)[^(l  -  x};(1  -  x/*)K g  ~  A;,)] 


[1*1 
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where 


K  =  exp 
a  ' 


K  p  =  exp 


f-K/'B  -  "X>  -  <"b  -  "a) 

TTT 


and 


Dj  —  ft/A“  exp  {- 


[19] 


[20] 


which  is  called  the  interface  diffusivity. 

In  order  to  maintain  steady  state  conditions  at  the  growing  interface  a  mass  balance 
must  be  satisfied.  In  order  to  extract  a  simple  model  for  the  present  we  will  assume  that 
interdiffusion  in  the  liquid  is  fast  enough  to  prevent  lateral  composition  variations  in  the 
first  liquid  monolayer  and  that  interdiffusion  in  the  solid  is  slow  enough  that  lateral 
exchanges  (ordering)  in  the  first  solid  monolayer  do  not  occur.  In  this  case  mass  balance 
requires  that 

JD  =  ?y(xB  ~  xfP 

JD  =  n^xB_xB^  t21l 


where  V  is  the  interface  velocity.  Equating  expressions  for  the  net  diffusive  fluxes  we 
obtain 


ft/,  Lx,/  L /,  aw  ^  _  V /  L  ft\ 

xB^1-xB^Ka  B^1  xB^Ka  “Vj^B  XB^ 


0(1  Lu<A/?_  xLfl„%(1  V_V  ,L_  0 

b(!-xb)K^  xB(l  xq)K^  -^-(Xg  Xg 


[22] 


where  Vg  =  D-/ A. 

One  may  choose  to  solve  these  equations  for  a  given  V  and  Xg  in  order  to  find  the 
compositions  of  the  a  and  0  sublattices.  This  would  yield  the  average  solid  composition 
and  the  long  range  order  parameter  at  the  interface.  For  this  purpose  we  require  the 
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redistribution  potentials  which  depend  on  temperature  and  the  various  compositions 
through  a  specific  thermodynamic  model  of  the  liquid  and  solid  phases.  To  determine  the 
temperature  of  the  interface  an  additional  equation  relating  the  solidification  velocity  to 
the  driving  force  for  solidification  is  required. 

4.2  The  Solidification  Velocity-Temperature  Relation 

We  will  assume  (28)  that  the  solidification  velocity  V  is  related  to  the  free  energy 
change  for  crystallization  of  AG  according  to 


v  =  vc 


AG  J 
TTT) 


[23] 


where  V^,  is  a  crystallization  rate  parameter.  Following  the  general  concept  that 
crystallization  is  easier  than  solute  redistribution,  Vq  is  assumed  greater  than  Vg.  In  the 
case  of  collision  limited  growth,  V^,  may  approach  the  speed  of  sound.  If  diffusive  jumps 
are  required  for  crystallization  Vq  can  be  approximated  by  the  ratio  of  the  liquid 
diffusion  coefficient  to  A,  the  jump  distance.  Thus,  given  thermodynamic  potentials  for 
the  phases  and  kinetic  parameters  and  V^,  Eqns.  [22]  and  [23]  give  the  order 
parameter,  liquid  composition,  and  temperature  at  the  interface  for  any  solid  composition 
and  solidification  velocity. 

5.  Choice  of  Thermodynamic  Potentials 

The  predictions  of  this  model  will  be  examined  for  phase  diagrams  having  a 

S  1 

congruently  melting  ordered  phase  at  Xg  =  Thus  the  order-disorder  temperature  of 
the  solid  phase  must  be  chosen  to  be  higher  than  the  congruent  melting  temperature.  The 
following  form  of  the  solid  free  energy,  G^,  will  also  permit  the  comparison  of  the 
predictions  when  this  order-disorder  transition  is  first-order  and  second-order: 

Gm  =  GAO-xB)  +  GBxE 
+  ni  {xl(! -x|)  +  ^,?2} 

+  n2r/4  [24] 
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+  \  R-T  |(xp  -  57?)  (n  (\p  -  i 7 f)  +  (Xg  +  2?/)  hi  (xjj  +  \i}) 

+  0  - ~hv)  fl>  (1  -  >4  - k'i)  +  U  _  xb  +  2*/) Cn  (l  _  xb  +  3r/)j 

s  s 

where  G^  and  Gg  are  the  free  energies  of  pure  solid  A  and  B  and  ft  j  and  ft9  arc 

constants.  If  ft9  =  0  and  ftj  <  0,  Eqn.  [24j  describes  a  Bragg-Williams  phase  with  near 

neighbor  interaction  ft  ^  and  ideal  entropy.  The  p-dependence  of  the  free  energy  curve  for 
S  1 

XB  =  2  IS  s^own  schematically  in  Figure  la.  The  transition  between  the  ordered  and 
disordered  phase  is  a  second  order  phase  transition;  i.e.,  the  equilibrium  order  parameter 
goes  to  zero  continuously  as  temperature  is  increased.  If  ft9  <  0  and  ftj  >  12ft9  the  fee 

s  1 

//-dependence  of  the  free  energy  curve  for  Xg  =  ^  is  shown  schematically  in  Figure  lb. 

The  transition  between  ordered  and  disordered  phase  is  a  first  order  phase  transformation; 
i.e.,  the  equilibrium  order  parameter  has  a  discontinuous  drop  to  zero. 

C 

For  the  choice  of  G  given  in  Eqn.  (24],  Eqns.  [9-12]  give 

"B  -  "A  =  °B  -  GA  +  ni<‘  -  2*B  - 

+  RT  (fu  (xg  -  j9j  -  fn  (I  -  x§  +  ji|)l 

#<B  -  (*a  =  °B  ”  GA  +  ni(1  "  2xB  +  1)  +  SR)’3 
+  HT  fto  (x|  +  U  -  ft.  (1  -  [25] 

M  +  "A>  =  °A  +  ni{(xI)2'J->2)  -»ty4 

+  ^RT(ft,[(l-xg)2-|,2]} 

+  /ig)  =  Gg  +  n,  {(1  -  Xg)2  -  ^r/2}  -  3rt2V4 

+  ^RT{to[(x|)2-|,2l)  [26] 

and 

(/*B  ~  V$f  =  GB  “  GA  +  Ql^  -  2xB  ~  V)  -  Sft^3 
(4-4)'  =  GB  -  GA  +  *Vl  -  2xB  +  9)  +  8Q2v3.  [27] 

For  the  liquid  phase  we  will  employ  a  regular  solution  model  with 
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;1A  =  Ga  +nL(x^)2  +  RTfr.(l-x^) 
/<[[  =  gJj  +  nL(l  -  Xq)2  +  RT  tn  xjj 


[28] 


where  G  ^  and  G ^  are  free  energies  of  pure  liquid  A  and  B.  The  redistribution  potential 


is  given  by 


0*B-^  =  GB"GA+ftL(1-2xB)' 


Further  we  will  assume  that 


GA-GA  =  GB-G!  =  SlTm-T> 


[29] 


[30] 


where  S  and  Tm  are  the  entropy  of  fusion  and  melting  point  of  the  pure  components 


(assumed  equal).  Eqn.  [30]  will  guarantee  a  phase  diagram  which  is  symmetric  about 

_  1 
XB  -  2‘ 

6-  Approximate  Result  for  Dependence  of  Order  Parameter  on  Velocity  at  a  Congruent 
Melting  Maximum 

A  simple  expression  can  be  obtained  from  Eqns.  [22]  for  x^  =  ^  which  will  provide 

some  insight  into  the  numerical  results  described  in  Section  7.  Both  equations  are 
S  1 

satisfied  if  Xg  =  7  and  if  the  order  parameter  satisfies 


+  sn2  v 

RT - 


=  tn 


1  +  v 


[31] 


for  t)  >  02.  Note  that  two  solutions,  77  =  0  and  77  #  0,  exist.  For  V  =  0  this  equation  gives 

C 

the  equilibrium  order  parameter  as  a  function  of  T  for  Xq  =  0.5  which  could  also  have 
been  found  using  Eqn.  [2].  As  a  kinetic  equation  for  solidification,  T  represents  the 
interface  temperature  which  can  only  be  determined  using  Eqn.  [23]  in  combination  with 
Eqn.  [31].  However,  if  we  consider  that  the  interface  temperature  may  be  approximated 


2For  i)  >  0,  Eqn.  [31]  is  modified  due  to  the  Aj  factors  and  the  theory  predicts  the  same 
behavior  for  ±  tj. 
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hv  the  melting  point  of  the  ordered  phase  T^,  then  Eqn.  [31],  after  rearrangement, 


becomes 


JljT?  +  812.-,  7/'" 


RT° 

m 


32 


n 


V  1 

This  approximation  predicts  that  77  goes  to  zero  as  y—  goes  to  (-  - — —  -  1 ).  For  the 

VD  2RT° 

case  wliere  120  =  0,  the  critical  temperature  for  the  order-disorder  transition  T  is  equal 
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1  V 

to  (-  T^).  Therefore  the  critical  value  of  y—  where  the  long  range  order  parameter  goes 


to  zero  is  given  by 


T 

/V  * 

D  T>  =  °  m 


[34] 


and  hence  the  critical  velocity  is  high  when  Tc  is  high  (i.e.,  when  the  ordering  is  strong). 
When  fi9  /  0,  the  critical  velocity  follows  the  same  trend  but  is  not  so  simply  estimated. 
7.  Numerical  Results 

We  will  examine  in  parallel  the  calculated  phase  diagram  and  results  of  the 
non-equilibrium  interface  theory  for  two  cases  with  thermodynamic  parameters  given  in 
Table  I. 


Case  I 
Case  II 


Table  I.  Choice  of  Thermodynamic  Parameters 
Tm(K)  S/R  ftL/R(K)  fyR(K)  0>/R(K) 

1000K  2  0  -3000  0 


1000K  2 


0  -2000  -566.2 


7.1  Equilibrium  Phase  Diagram 

Figures  3a  and  3b  show  the  equilibrium  order  parameter,  j;  ,  of  the  solid  as  a 
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c 

function  of  T  for  the  two  cases  for  Xg  =  0.5,  obtained  using  Eqn.  [2].  In  I,  the 

n. 

order-disorder  (OD)  temperature  is  1500  K  (-  In  II,  the  OD  temperature  is  1583  K 
and  is  obtained  by  finding  the  temperature  where  Gg  for  r?  =  0  is  equal  to  Gg  for  r/  =  r)Q. 
The  dashed  curve  represents  states  where  the  ordered  phase  is  metastable  whereas  the 
dotted  curve  represents  states  where  the  ordered  phase  is  unstable.  Cases  I  and  II  are 
respectively  second  and  first  order  transitions. 

Figures  4a  and  4b  show  half  of  the  phase  diagrams  for  the  solid  phases  only.  In 
Figure  4a  no  two-phase  field  exists  between  the  ordered  and  disordered  phase, 
characteristic  of  the  second  order  transition.  The  curve  is  obtained  from  Eqn.  [2]  by 
letting  ;/  -♦  0.  In  Figure  4b  a  two-phase  field  exists.  The  boundaries  are  obtained  by 
equating  the  values  of  ^(/i“  +  for  T)  =  0  and  V  =  obtained  from  Eqn.  [2]  and 
similarly  the  values  of  ^(/*g  +  /*g)- 

Figures  5a  and  5b  show  the  equilibrium  diagrams  using  Eqns.  [3]  in  combination 
with  Eqn.  [2]  to  generate  the  liquid-solid  equilibrium.  In  both  cases  the  melting  point  of 
the  composition  .\g  =  0,5  is  below  the  OD  temperature.  The  diagrams  contain  "peritectic 
reactions,"  L  +  O  -*  D;  however,  in  Figure  5a  no  discontinuity  in  slope  occurs  at  the 
liquidus  at  the  "peritectic  temperature"  due  to  the  second  order  nature  of  the  OD 
reaction.  For  case  II  a  metastable  liquidus  and  solidus  are  shown  dashed  for  the 
disordered  phase.  For  case  I  such  curves  would  represent  unstable  equilibrium  and  are 
not  shown. 

Xs 

7.2  Non-Equilibrium  Interface  Conditions  for  Coneruentlv  Melting  Ordered  Phase 

In  general  there  are  two  solutions  to  Eqns.  [22]  and  [23]  for  the  non-equilibrium 
interface  conditions:  one  for  i?  =  0  and  one  for  tj  i  0.  The  rj  i  0  solution  only  exists  below 
some  critical  value  of  V/Vp.  Different  parts  of  each  solution  are  stable  or  unstable  in  the 
sense  that  small  fluctuations  in  the  order  parameter  will  lead  to  large  changes  in  the 
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interface  conditions.  Unstable  portions  will  be  shown  dotted  in  the  following  figures.  The 


details  of  this  stability  analysis  will  be  deferred  until  the  discussion. 

For  4  =  0.5  Eqns.  [22]  and  [23]  are  satisfied  for  =  0.5  and  values  of  q  and  T  as 

VC  :j 

a  function  of  A—-  shown  in  Figures  6  and  7.  For  these  results  t was  chosen  to  be  10  . 
v  D  D 


In  Figure  6a  for  case  I,  the  order  parameter  goes  to  zero  continuously  as  a  critical  value  of 

X—  is  reached.  In  Figure  6b  for  case  II  the  order  parameter  exhibits  a  turning  or  limit 
VD 

point  at  a  critical  value  of  X—.  If  the  velocity  is  increased  above  the  limit  point,  the 

VD 


order  parameter  must  jump  discontinuously  to  the  7?  =  0  solution  in  case  II.  In  Figure  6b 

the  branch  below  the  turning  point  is  dotted  and  represents  unstable  solutions.  At  V  =  0 

this  unstable  branch  corresponds  to  the  thermodynamically  unstable  equilibrium  between 

the  liquid  and  solid  (see  T  =  in  Figure  lb). 

In  Figure  7,  the  interface  temperatures  for  the  q  —  0,  and  q  t  0  solutions  are  shown. 

Again  the  dotted  portions  are  unstable.  For  case  I,  the  q  —  0  solution  is  only  stable  above 

the  critical  velocity.  For  case  II  the  q  —  0  solution  is  always  stable  and  corresponds  to 

the  solidification  of  the  metastable  disordered  phase. 

Figure  8  shows  cross  plots  of  q  and  T.  Points  along  these  curves  are  parameterized 

by  X— .  Also  included  in  the  figure  are  the  equilibrium  order  parameter,  q  ,  versus 
D  e 


temperature  as  shown  originally  in  Figure  3. 

7.3  Simultaneous  Disorder  and  Solute  Trapping 

For  compositions  away  from  the  congruent  maximum  of  the  ordered  phase,  the 
possibility  of  non-equilibrium  incorporation  of  solute  (solute  trapping)  as  well  as  disorder 
trapping  is  possible.  Figures  9  and  10  show  the  order  parameter  and  liquid  composition  at 

c 

the  interface  as  a  function  of  V/V^  for  xR  -  0  .48  for  the  two  cases.  The  behavior  of  the 

c 

order  parameter  is  quite  similar  to  that  for  x^  =  0.5.  Figure  10  is  shown  with  a 
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logarithmic  scale  for  velocity  to  permit  comparison  to  previous  theory  (1-2)  on  solute 
trapping  of  disordered  solids.  In  both  cases  I  and  II  the  composition  of  the  liquid 
approaches  the  composition  of  the  solid  at  high  velocity. 

Near  the  velocity  where  rj  goes  to  zero  in  case  I  or  where  the  turning  point  occurs  in 
case  II  in  Figure  9,  the  liquid  compositions  shown  in  Figure  10  have  abrupt  changes.  In 
case  I  the  lower  curve  for  the  r;  ^  0  solution  merges  with  that  for  the  q  =  0  solution. 

Below  this  merge  point  the  q  =  0  solution  is  unstable  and  is  shown  dotted.  In  case  II,  the 
liquid  composition  curve  for  the  q  -  0  solution  is  stable  for  all  velocities.  The  liquid 
composition  for  the  q  /  0  solution  exhibits  a  turning  point  at  the  same  velocity  where  q 
has  a  turning  point  in  Figure  9b.  In  this  case  the  lower  branch  is  stable  and  the  upper 
branch  is  unstable.  Thus  for  case  II  the  theory  predicts  that  with  increasing  velocity  a 
jump  in  liquid  composition  at  the  interface  occurs  when  ??goes  to  zero.  As  shown,  the  size 
of  the  jump  is  small  but  it  of  course  depends  on  the  thermodynamic  potentials  used. 

It  is  also  interesting  to  note  that  in  case  I  q  goes  to  zero  at  velocities  where  solute 
redistribution  remains  significant  while  in  case  II  complete  disorder  trapping  occurs  at  a 
velocity  where  solute  trapping  is  practically  complete.  Numerical  results  for  other  values 
of  Qj  and  not  described  here  suggest  that  this  is  not  a  general  trend  for  solid  phases 
with  first  and  second  order  transitions. 

8.  Discussion 

The  validity  of  the  assumptions  regarding  complete  mixing  in  the  first  liquid 
monolayer  and  the  absence  of  lateral  atomic  rearrangements  in  the  first  solid  monolayer 
are  ultimately  tied  to  the  values  of  the  interdiffusion  coefficients  in  the  solid  (Dg), 
interface  (Dj ),  and  liquid  (DjJ.  In  general  it  is  reasonable  to  assume  that  Dg  <  Dj  <  Dj . 
The  present  theory!  best  applies.iwhen  these  differ  by  orders  of  magnitude.  However,  even 
if  Dj  approaches  Dj ,  the  assumption  of  complete  mixing  in  the  first  liquid  monolayer  may 
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only  become  suspect  if  the  atom  positions  of  each  sublattice  stays  fixed  laterally  on  the 
liquid  solid  interface  during  growth.  This  could  happen  in  some  crystallographic 
orientations  of  some  crystal  structures;  e.g.,  growth  of  a  chemically  layered  structure 
edgewise  into  the  melt.  In  this  case  one  would  expect  disorder  trapping  at  lower 
solidification  rates  than  predicted  by  this  theory.  Relaxation  of  the  assumption  regarding 
no  ordering  in  the  first  solid  monolayer  would  increase  the  state  of  order  at  any  velocity 
over  that  predicted  by  the  theory.  This  ordering  of  the  surface  layer  could  be  added  to 
that  which  occurs  during  solid  state  cooling  to  predict  the  state  of  order  in  a  final  (cold) 
solid.  The  theory  also  has  assumed  that  or  that  the  interdiffusion  between  the 

liquid  and  solid  do  not  depend  on  the  sublattice.  This  condition  can  easily  be  relaxed. 

One  result  of  such  a  change  would  be  that  for  solidification  of  a  composition  at  a 

L  S 

congruent  melting  maximum,  x^  would  not  equal  x^  except  at  V  =  0  and  V  =  ®. 

The  results  of  the  dependence  of  tj  on  X—  may  be  viewed  in  the  context  of 

VD 


bifurcation  theory.  Bifurcation  theory  describes  the  change  in  stability  of  solutions  near 

intersection  (bifurcation)  points  and  turning  points.  The  two  solutions  r)  =  0  and  r;  f  0 

V  **  1  \; 

always  intersect  at  -n—  ~ - —  -  1,  which  occurs  at  positive  and  negative  values  of  r- — . 

n  2RT°  v  D 

for  case  I  and  II  respectively.  At  this  point  the  stability  of  the  77  =  0  solution  must 

change.  By  replacing  Eqn.  [31],  which  applies  only  for  the  steady  state,  with  the  time 

dependent  form  of  the  flux  balance,  erne.  obfauin-S  fo  r 


dt  =~?{2VT  +  0  +  9)  exp 


D 


n j  77  +  -1 


RT 


[35] 


m 


,  ,  ,  Jor  ^>0’ 

frvtn  uhich 

^  the  stability  of  the  various  solutions  can  be  inspected.  If  a  positive  perturbation  of  r? 
makes  the  RHS  of  Eqn.  [35]  positive  then  ^  <  0  and  the  perturbation  decays  and  that 
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portion  of  the  solution  is  stable.  Using  this  method  the  77  =  0  solution  can  be  shown  to  be 
V  ^  1 

stable  only  if  rj—  >  -1 - - — .  Similarly  the  stability  of  the  77  f  0  solutions  can  be 

VD  2RT° 

inspected.  In  case  1,  the  7/  /  0  solution  is  always  stable  while  the  77 1  0  solution  for  case  II 


must  change  from  stable  to  unstable  at  the  turning  point.  The  stable  branch  has  negative 


slope  in  the  77  versus  y—  plot.  3r 

-  r  • .  im 

The  slope  of  the  curves^in  Figure  7,  are  worthy  of  some  discussion.  The  slope  is 
called  the  interface  kinetic  coefficient  //  and  relates  a  simple  measure  of  the  departure 


from  equilibrium  required  to  cause  growth  at  various  rates.  It  can  be  seen  that  the  slope 


is  small  near  V  =  0  and  whenever  77  =  0.  At  velocities  where  77  is  changing,  the  slope  is 

d  T 

much  steeper;  i.e.,  the  kinetics  become  more  "sluggish."  A  complex  expression  for  £jy  can 
be  found  for  Xg  =  ^  by  differentiation  of  Eq.  [23]  and  [24]. 

RTm  1 

Near  V  =  0  or  when  77  =  0,  the  kinetic  coefficient  is  — ^^{y— )  where  AS  is  the  plicable 

^  v. .  ~  : 

entropy  of  fusion.  When  77  is  changed  rapidly  the  kinetic  coefficient  depends^on  both  VC 
rvi  io  tMkI  Vg.  Thus  the  growth  kinetics  are  controlled  by  either  when  the  solid  is 

disordered  or  when  the  solid  is  ordered  at  small  V.  At  intermediate  V  the  kinetics  are 


controlled  by  a  combination  of  and  Vp  when  the  solid  is  ordered. 

The  present  theory  can  provide  some  guidance  in  selecting  intermetalliccompounds 
which  respond  to  rapid  solidification  to  produce  phases  with  non-equilibrium  composition 
or  state  of  order.  The  former  may  be  useful  to  provide  the  supersaturation  necessary  to 
form  a  controlled  scale  of  precipitation  by  subsequent  heat  treatment.  Non-equilibrium 
order  produced  by  solidification  will  often  be  followed  by  rapid  ordering  during  solid  state 
cooling.  The  fine  antiphase  domains  structures  produced  by  this  process  can  provide 
unusual  starting  materials  for  subsequent  heat  treatment  (30). 
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It  is  important  to  note  that  many  complex  ordered  phases  such  as  a  phases  bear  no 
resemblance  to  simple  disordered  structures.  The  elimination  of  such  phases  by  rapid 
solidification  and  replacement  by  a  simpler  disordered  phase  is  not  an  example  where  this 
theory  is  applicable.  Classical  concepts  of  nucleation  and  growth  competition  must  be 
applied  to  these  situations. 

Quantitative  predictions  by  this  theory  require  good  thermodynamic  models  for  the 
phases  of  interest  as  well  as  values  for  the  kinetic  parameters  V ^  and  Vp.  Unfortunately, 
measurements  of  equilibrium  order  parameter  are  often  not  available,  especially  for 
materials  of  technological  interest  with  high  melting  points.  Thus  the  data  base  for 
thermodynamic  modelling  is  poor.  Additionally  functional  forms  of  the  thermodynamic 
potentials  used  for  ordered  phase  in  many  phase  diagram  evaluations  have  improper 
functional  form  near  77  =  0.  This  will  severely  impact  the  predictions  of  this  theory. 

Values  of  Vc  and  can  only  be  crudely  estimated.  Experiments  on  combined  solute 
trapping  and  disorder  trapping  using  pulsed  laser  annealing  with  the  transient 
conductance  method  to  measure  interface  rates  may  provide  insight  in  these  parameters. 
An  analysis  of  the  rate  of  ordering  during  solid  state  cooling  must  be  included  to  obtain 
the  values  produced  by  the  interface  itself.  Two  features  are  amenable  to  measurement. 
Extremely  fine  antiphase  domains  should  be  present  if  the  critical  velocity  has  .been 
exceeded.  Also  the  abrupt  behavior  of  the  composition  shown  in  Figure  10  when  t)  goes  to 
zero  may  produce  unusual  microstructural  features. 
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Figure  Captions 

Figure  1.  Schematic  representation  of  the  origin  of  the  thermodynamic  constraint  on  7? 
during  non-equilibrium  solidification  of  a  congruently  melting  ordered  phase.  The 

K'«-/ 

ordered  solid  phase  differ/ from  the  disordered  solid  phase  by  (a)  a  second  order  transition 

O  Y 

■asd  (b)  a  first  order  phase  transition.  Top,  solid  free  energy  versus  77  and  liquid  free 
energy  varies  temperatures  Tp  T2,  T^.  Tj  is  the  equilibrium  melting  point  of  the 
ordered  phase,  T9  is  an  arbitrary  temperature  for  (a)  and  the  metastable  melting  point  of 
the  disordered  phase  for  (b),  and  T^  is  the  unstable  melting  point  of  the  ordered  phase. 
The  condition,  Gg  -  G^  <  0,  defines  a  range  of  77  as  a  function  of  solidification 
temperature  where  solidification  is  possible  at  the  bottom. 

Figure  2.  Reaction  coordinate  diagram  for  interface  redistribution  reaction.  Initial  states: 
B  on  a  or  p  site  in  solid,  A  in  liquid.  Final  states:  A  on  a  or  0  site  in  solid,  B  in  liquid. 
Figure  3.  Calculated  equilibrium  long  range  order  parameter  in  the  solid,  77,  as  a  function 

C 

of  temperature  for  Xg  =  0.5  (a)  Case  I,  order-disorder  temperature  is  1500  K.  (b)  Case 
II,  order-disorder  temperature  is  1583  K,  dashed  curve  -  metastable  ordered  state,  dotted 
curve  unstable  ordered  state. 

Figure  4.  Calculated  phase  diagram  for  solid  phase  only  for  (a)  Case  I  and  (b)  Case  II. 
(Diagram  is  symmetric  about  Xg  =  0.5). 

Figure  5.  Calculated  equilibrium  phase  diagram  for  (a)  Case  I  and  (b)  Case  II.  In  (b)  the 
dashed  curve  corresponds  to  the  metastable  liquidus  and  solidus  for  the  disordered  phase. 
Figure  6.  Long  range  order  parameter  of  the  solid  at  the  interface,  77,  as  a  function  of 

S 

dimensionless  growth  velocity,  V/Vg  for  the  congruently  melting  compound  at  Xg  =  0.5. 
(a)  Case  I,  (b)  Case  II.  Capital  letters  designate  the  common  interface  conditions  in 
Figures  6,  7,  and  8. 

Figure  7.  Interface  temperature  T  as  a  function  of  dimensionless  growth  velocity,  X—  for 

V  rv 


23 


the  eongruently  melting  compound  at  x^  —  0.5.  (a)  Case  I,  (b)  Case  II.  The  dashed 
curves  correspond  to  the  //  =  0  solution. 

Figure  S.  Cross  plot  of  order  parameter  and  interface  temperature  parameterized  by  X— 

N  I) 

for  (a)  Case  I  and  (b)  Case  II.  Also  shown  are  the  equilibrium  order  parameter  at  the 
various  temperatures  taken  from  Figure  3. 

Figure  9.  Long  range  order  parameter  of  the  solid  at  the  interface,  r],  as  a  function  of 

C 

dimensionless  growth  velocity,  V/V^  for  (a)  Case  I  with  =  .48,  and  (b)  Case  II  with 

S 

Xrj  =  0.48.  Capital  letters  designate  common  interface  conditions  in  Figures  9  and  10. 
Figure  10.  Liquid  composition  at  the  interface  x^,  as  a  function  of  dimensionless  growth 
velocity  V/V^  same  cases  as  Figure  9. 
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